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Анотація. У роботі розглядається інтелектуальна система комп'ютерної математики Іпрагбоїуег, яка 
призначена для автоматичного дослідження та розв'язування основних класів задач обчислювальної матема- 
тики на багатоядерних комп'ютерах з графічними прискорювачами. Окреслено проблеми достовірності ре- 
зультатів розв'язання задач із наближеними даними. Проаналізовано особливості використання існуючих сис- 
тем комп'ютерної математики, виявлено їх недоліки. Описуються функціональні можливості Шраг50Їуег, деякі 
інноваційні підходи щодо реалізації ефективного розв'язування задач на гібридній архітектурі. Наведено прик- 
лади застосування Шраг8оЇуег у математичному моделюванні процесів у різних предметних областях. У наш 
час, у багатьох предметних областях (ядерна енергетика, механіка, хімія, молекулярна біологія, медицина 
тощо) постійно виникають нові, більш складні об'єкти та явища, що підлягають математичному дослідженню 
на комп'ютері. Це спонукає до розробки нових чисельних методів та технологій математичного моделювання, 
а також до створення більш потужних комп'ютерів для їх реалізації. Із появою та постійним розвиненням 
суперкомп'ютерів різної архітектури, дуже великої актуальності набувають проблеми їх ефективного викорис- 
тання, розширення кола задач, які необхідно розв'язувати, забезпечення достовірності комп'ютерних результа- 
тів та підвищення рівня інтелектуальної інформаційної підтримки користувачів - фахівців у різних предметних 
областях. Питанням вирішення цих проблем приділяється на сьогодні особлива увага багатьох спеціалістів у 
галузях інформаційних технологій та паралельного програмування. Найвидатніші в світі вчені в галузі комп'ю- 
терних технологій вирішення проблем ефективного використання сучасних суперкомп'ютерів вбачають у 
створенні алгоритмічно-програмного забезпечення, яке легко адаптується до різних комп'ютерних архітектур з 
різними видами пам'яті та співпроцесорами, підтримує ефективний паралелізм на мільйонах ядер тощо. Крім 
того, підвищення ефективності реалізації високопродуктивних обчислень на сучасних суперкомп'ютерах 
передбачається шляхом їх інтелектуалізації, передаючи комп'ютеру для самостійного виконання значну части- 
ну функцій (символьні мови для комп'ютерної постановки задачі, дослідження властивостей математичних 
моделей, візуалізація та аналіз результатів розв'язування задач тощо). Індустрія розробки та використання 
інтелектуальних комп'ютерних технологій є одним із основних напрямків розвитку науки та технологій у 
сучасному суспільстві. 

Ключові слова: інтелектуальна система комп'ютерної математики; математичне моделювання; набли- 
жені дані; засоби штучного інтелекту. 
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Вступ 

Одним з перспективних шляхів вико- 
ристання високопродуктивних комп'юте- 
рів при математичному моделюванні про- 
цесів з великими обсягами даних в науці та 
інженерії на сьогоднішній день є викорис- 
тання комп'ютерів гібридної архітектури, 
яке поєднує МІМЮ-архітектуру (багато- 
процесорні, багатоядерні СР) та комп'ю- 
тери із 51МО-архітектурою графічних при- 
скорювачів (СРО), - гібридні комп'ютери. 
Проте, при використанні цих комп'ютерів 
у математичному моделюванні виникають 
проблеми ефективного використання двох 
різних комп'ютерних архітектур. 

Як відмічають багато спеціалістів, 
досить великий набір існуючого програм- 
ного забезпечення для математичного мо- 
делювання, як правило не враховує похиб- 
ки вихідних даних розв'язуваної задачі (|1- 
11). Проте, математичні моделі фізичних 
процесів, які описуються диференціальни- 
ми та інтегральними рівняннями або варіа- 
ційними задачами, можуть мати природні 
заокруглення в результаті вимірювань, спо- 
стережень, припущень, гіпотез 1 т. п. У цьо- 
му випадку вихідні дані є деяким набли- 
женням до точних даних задачі. 

Надалі, при дискретизації матема- 
тичної моделі, ці похибки трансформують- 
ся у похибки коефіцієнтів рівнянь, які під- 
лягають розв'язуванню. Крім того, вихідні 
дані математичних моделей можуть бути 
задані точно в числовому вигляді або 
представлені математичними формулами, 
але заокруглення при введенні чисел в 
комп'ютер та обчислення формул також 
призводить до комп'ютерних моделей з 
наближеними даними |З, 10, 12, 131. 

Проблеми достовірності комп'ютер- 
них результатів розв'язування задач з на- 


ближеними даними можуть бути вирішені 
в тому випадку, якщо в комп'ютерному 
середовищі дослідити математичні влас- 
тивості комп'ютерної моделі задачі та роз- 
в'язати алгоритмом у відповідності до цих 
властивостей з аналізом достовірності 
отриманих результатів, враховуючи архі- 
тектурні та технологічні особливості 
комп'ютера ПО). 

Самостійне вирішення вказаних 
проблем математичного моделювання на 
гібридних комп'ютерах вимагає від 
користувачів-спеціалістів з різних пред- 
метних областей відповідних знань в об- 
ласті обчислювальної математики, архі- 
тектурних особливостей комп'ютерів, на- 
виків використання різних технологій па- 
ралельного програмування 1 т. д. 

У Інституті кібернетики 
імені В. М. Глушкова НАН України ство- 
рено інтелектуальну систему комп'ютер- 
ної математики (ІСКМ) Праг5оїуег для ав- 
томатичного виконання на гібридних ком- 
п'ютерах дослідження та розв'язування (з 
гарантією достовірності результатів) задач 
обчислювальної математики з наближени- 
ми даними, які виникають при математич- 
ному моделюванні складних процесів та 
явищ у наукових та інженерних до- 
слідженнях, а саме: системи лінійних ал- 
гебраїчних рівнянь, алгебраїчна проблема 
власних значень, системи нелінійних рів- 
нянь, задачі Коші для систем звичайних 
диференціальних рівнянь з наближеними 
даними. При цьому автоматично викону- 
ється розпаралелення обчислень та даних 
при ефективному використанні обчислю- 
вальних ресурсів гібридного комп'ютера. 
Інтелектуальний інтерфейс надає можли- 
вість вводити дані задачі на мові предмет- 
ної області та отримувати результати роз- 
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в'язування з детальним поясненням про- 
цесу дослідження задачі. Таким чином, 
Іпраг5оімег забезпечує роботу користувача 
на комп'ютері зі складною гібридною ар- 
хітектурою як на комп'ютері однопроце- 
сорної архітектури. 

Тут ми будемо розглядати універ- 
сальні функціональні можливості Іпраг- 
Зоїуег на класі задач - системи лінійних 
алгебраїчних рівнянь (СЛАР). 

Аналіз проблем використання існу- 
ючих систем комп'ютерної математики 

У наш час велика увага приділяється 
розробці та використанню систем комп'ю- 
терної математики (СКМ). Це поняття 
включає сукупність як теоретичних та 
технологічних засобів, так 1 сучасних 
програмних й апаратних засобів, що нада- 
ють можливість: формулювати завдання 
за допомогою спеціальної символьної мо- 
ви в термінах математики; автоматизувати 
виконання як чисельних, так і аналітичних 
обчислень; розробляти складні обчислю- 
вальні алгоритми для розв'язування задач; 
застосовувати засоби візуалізації процесів 
та даних тощо | 141. 

Найбільшу популярність мають такі 
СКМ: МаїСаАР, Маріе, Майетайса, 
Мацар. Ці програмні засоби найчастіше 
використовуються для розв'язання на- 
вчальних, інженерних, науково-дослідних 
задач у різних галузях природничих наук. 
Серед недоліків цих СКМ користувачі від- 
мічають, що реалізація чисельних алго- 
ритмів (наприклад, розв'язування неліній- 
них або диференціальних рівнянь) в них 
базується на точно заданих даних |1, 15). 
Тому є випадки отримання неправильних 
результатів розв'язування, не надається 
інформація, яким методом розв'язувалася 
задача. Документація для користувачів до- 
сить велика (наприклад, обсяг документа- 
ції МаЧаб досягає майже 5 тис. сторінок, 
що робить її важко доступною для огля- 
ду). Специфічні редактори коду програм 
також ускладнюють роботу з СКМ. 

За десятиріччя використання ком- 
п'ютерів для розв'язування широкого кола 
різноманітних прикладних задач створено 
велику кількість проблемно орієнтованих 
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програмних пакетів, наприклад, АМ5У5, 
МАЗТКАХМ та інші. В переважній більшос- 
ті ці програмні засоби призначено для од- 
нопроцесорних комп'ютерів. Для ефек- 
тивного використання на суперкомп'юте- 
рах, як відмічають спеціалісти |Ї, 2, 5-9|, 
вони потребують суттєвої інтелектуаліза- 
ції - створення засобів автоматичного на- 
лаштування на архітектуру комп'ютера та 
автоматичної побудови математичних мо- 
делей і відповідних паралельних алгорит- 
мів їх дослідження, а також розробки ін- 
терфейсів для спілкування з користувача- 
ми на мові предметної області. 

Таким чином, створення ІСКМ з 
функціями автоматичного дослідження та 
розв'язування задач обчислювальної мате- 
матики та гарантією достовірності резуль- 
татів є актуальною роботою. 

Гібридні алгоритми для досліджен- 
ня та розв'язування СЛАР в Праг5оіуег 

Для дослідження й розв'язування лі- 
нійних систем на гібридних комп'ютерах 
створено низку нових паралельних алго- 
ритмів відомих прямих методів для щіль- 
них та розріджених матриць різної струк- 
тури на основі Щ/, / та 5УД-розвинень 
матриць (10, 16, 17). Також представлено 
нові гібридні алгоритми та програми іте- 
раційних методів розв'язування СЛАР з 
розрідженими матрицями  нерегулярної 
структури та наближеними даними на 
основі | методу верхньої релаксації, 
поперемінно трикутного методу тощо 
15, 19). 

Із метою підвищення ефективності 
розв'язування задач лінійної алгебри для 
розріджених матриць нерегулярної струк- 
тури великих розмірів використано ідею 
попереднього їх зведення за допомогою 
методу паралельних перерізів до блочно 
діагональної матриці з обрамленням |171. 
Таке представлення розрідженої матриці 
дає можливість зменшувати обсяг даних 
задач та більш ефективно виконувати роз- 
паралелення обчислень як на СР), так і на 
СРО. 

Кожен гібридний алгоритм здійснює 
дослідження його відповідності математич- 
ним властивостям СЛАР з наближеними 
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даними, які виявляються в комп'ютері, та 
розв'язування на ефективній топології ком- 
п'ютера з оптимальної кількості процесо- 
рів СРО та відповідних СРО. Процесори 
СРО, здебільшого, використовуються для 
багатопоточного виконання  матрично- 
матричних та матрично-векторних опера- 
цій при асинхронному копіюванні необхід- 
них даних великих обсягів між СРО та 
СРО. 

У разі необхідності уточнення ре- 
зультатів розв'язування, процес обчислень 
автоматично виконується з підвищеною 
комп'ютерною розрядністю. По закінчен- 
ню обчислювального процесу видається 
розв'язок задачі з оцінками точності або 
пояснення причини неможливості отрима- 
ти розв'язок. 

Одна із важливих проблем, яка вини- 
кає при розв'язанні практичних задач з 
розрідженими матрицями апріорі не виз- 
наченої структури, полягає у виборі ефек- 
тивного алгоритму. У ІСКМ для задач з 
розрідженими матрицями великих обсягів 
реалізовано за допомогою методів штуч- 
них нейронних мереж автоматичне розпіз- 
навання структури матриць та вибір най- 
ефективнішого гібридного алгоритму роз- 
в'язування за такою схемою |20): 

е візуалізація матриці на основі класифі- 
кації об'єктів; 

е- визначення типу («портрету») матриці; 

е визначення необхідного алгоритму роз- 
в'язування для типу матриці; 

е дослідження та розв'язання задачі від- 
повідним гібридним алгоритмом. 

Такий підхід до створення гібридних 
алгоритмів забезпечує його масштабова- 
ність, ефективне використання обчислю- 
вальних ресурсів та часу розв'язування, 
гарантує достовірність результатів розв'я- 
зування. 

Комп'ютерне дослідження матема- 
тичних властивостей СЛАР з наближе- 
ними даними 

Коротко опишемо методологію ком- 
п'ютерного дослідження математичних 
властивостей СЛАР з наближеними дани- 
ми в Шраг5оіуєг |З, 10). 


Розв'язування сумісних СЛАР 

Ах аб 0), 
де, в загальному випадку, А - прямокутна 
матриця розміру тхп; Б - матриця правої 
частини розміру тха (т-мірний вектор), 
зводиться до знаходження такого розв'яз- 
ку х (матриці розміру пхи, п-мірного век- 
тора), щоб рівняння (1) перетворювалося в 
тотожність. Відомо, якщо матриця А сис- 
теми квадратна, невироджена, (тобто Її 
визначник |А| є еК(А) 0), то розв'язок 
СЛАР (1) існує й єдиний. 

Несумісні системи класичного роз- 
в'язку не мають, проте, в цьому випадку, 
можна знайти такий розв'язок х, який мі- 
німізує евклідову норму | Ах - Ь| . Вектор х 
називається розв'язком за методом най- 
менших квадратів або узагальненим роз- 
в'язком СЛАР. У загальному випадку іс- 
нує нескінченна множина узагальнених 
розв'язків. Розв'язок, який має найменшу 


норму ТР 


, Єдиний і називається нормаль- 


ним узагальненим розв'язком або нор- 
мальним псевдорозв'язком. Іноді ставить- 
ся задача визначення стійкої проєкції на 
підпростір, утворений сингулярними век- 
торами матриці А. 

Як вже зазначалося, при розв'язу- 
ванні прикладних задач рідко виникають 
СЛАР з точними даними 

Ахар. (2) 
Найбільш типова є постановка задачі (1) 
разом із заданням відповідних похибок у 
вихідних даних: 


іл -АЇ-АЦє єдЇЙ| 


І» -5- МАВ ва?) 

При цьому передбачається, що структура 
матриці вихідної задачі (2) та збуреної за- 
дачі (1), (3) не змінюється, тобто, якщо ви- 
хідна матриця є симетрична, то й збурена 
залишається також симетричною, якщо 
вихідна - стрічкова, то й збурена - 
стрічкова, 1 т. д. 

Тут ми розглядаємо випадок, коли 
А -- квадратна матриця розміру пхи; 
Ь -- вектор правої частини СЛАР розміру 
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п; х - розв'язок (вектор розміру п); єд, 
Єр- максимальні відносні похибки еле- 


ментів матриці 1 правої частини. 

Розв'язуючи СЛАР із наближеними 
даними, необхідно розглядати цілий клас 
систем рівнянь (1), (3), що має досить ши- 
року множину формально допустимих 
розв'язків. Тому, розв'язавши задачу (1), 
необхідно оцінити збурення розв'язку, в 
залежності від збурення вихідних даних 
(3). Не завжди схожість елементів матриць 
А і А та правих частин Б і Б забезпечує 
достатню схожість розв'язків. Наприклад, 
при деякому збуренні в межах точності за- 
дання елементів матриці і/або правої час- 
тини несумісної, точно заданої системи 
(2), отримана в комп'ютері збурена систе- 
ма (1), (3) може виявитися сумісною й 
навпаки - сумісна СЛАР може перетвори- 
тися в несумісну | 10). 

У першу чергу, комп'ютерне до- 
слідження СЛАР з наближеними даними 
визначає існування та єдиність класичного 
розв'язку, сумісність системи. Крім того, 
досліджується коректність постановки за- 
дачі (1), (3), виродженість і/або додатно 
означеність матриці системи в межах по- 
хибки вихідних даних дослідження СЛАР 
з наближеними даними. Невиродженою в 
межах точності задання вихідних даних 
вважається матриця, яка не може стати ви- 
родженою в області ДА зміни числових 
значень її елементів з умови (3). 

Машинно невиродженою вважається 
матриця, яка не може стати виродженою 
при зміні числових значень її елементів у 
межах машинної точності. 

Теоретичним критерієм коректності 
СЛАР (1), (3) з квадратною матрицею та з 
наближеними вихідними даними є вико- 
нання умов: 


фекА)ж о, | |ЛАЦИ «1 


для довільного збурення ДА в межах за 
(3), які гарантують існування, єдиність та 
стійкість класичного розв'язку задач в об- 
ласті змін вихідних даних задачі (1), (3). 
Система з прямокутною матрицею 
апріорі вважається некоректною. Для сис- 
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тем з прямокутними та виродженими мат- 
рицями обчислюється нормальний уза- 
гальнений розв'язок, який, як вказувалося 
вище, у випадку сумісної системи співпа- 
дає із класичним розв'язком з найменшою 
нормою. 

Комп'ютерне дослідження корект- 
ності зводиться до перевірки двох спів- 
відношень: 


і0чугі0,  (ЦААЦАЇ)ЛСА) «1, (4) 
де П(А) - число обумовленості матриці 


системи, у - А). Перша умова, яка ви- 
конується в арифметиці із плаваючою ко- 
мою, означає, що матриця невироджена в 
межах машинної точності (машинно неви- 
роджена), а друга - що вона невироджена 
в межах точності задання вихідних даних. 

При виконанні умов (4) розв'язок 
комп'ютерної задачі існує, єдиний 1 стій- 
кий. Таку комп'ютерну задачу слід роз- 
глядати як коректно поставлену в межах 
точності введення вихідних даних. У ін- 
шому випадку, матриця системи може ви- 
явитися матрицею неповного рангу й ком- 
п'ютерну задачу (1), (3) слід розглядати як 
некоректно поставлену. 

Таким чином, основним критерієм 
для визначення властивостей СЛАР з не 
виродженими матрицями, які впливають 
на достовірність розв'язку задачі, є число 
обумовленості И(А) матриці системи. 

Для обчислення числа обумовленості 
невироджених матриць використовується 


обернена матриця А" 
пл) -ДАЦАГ'Ї. (5) 


Якщо значення Л(А) у формулі (35) неве- 
лике, то матриця системи називається доб- 
ре обумовленою, в іншому випадку - по- 
гано обумовленою. Проте, в комп'ютер- 
них обчисленнях критерій доброї або по- 
ганої обумовленості матриці є величиною 
відносною, пов'язаною з математичними 
можливостями конкретного комп'ютера, 
похибками даних задачі тощо. Одна й та 
сама матриця може класифікуватися для 
однієї довжини мантиси машинного слова 
як машинно добре обумовлена, а для ін- 
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шої - як машинно погано обумовлена. 

Погана «машинна обумовленість» 
тісно пов'язана із математичними можли- 
востями конкретного комп'ютера. Збіль- 
шуючи довжину машинного слова, можна 
одержати комп'ютерний розв'язок, який 
буде досить близьким до математичного 
розв'язку задачі. 

Обчислення оберненої матриці по- 
требує великих витрат часу та пам'яті 
комп'ютера, тому в Іпраг5оімег для до- 
слідження достовірності отримуваних роз- 
в'язків СЛАР з невиродженими матриця- 
ми використовується оцінка числа обу- 
мовленості, позначимо як сопад, яка об- 
числюється значно простіше (10). 

Якщо значення сопад в комп'ютері 
із плаваючою комою задовольняє умові: 


1.0--1.0/сопад «1.0, (6), 


то матриця в Іпраг5оіуег вважається ви- 
родженою в межах машинної точності й 
розв'язування СЛАР автоматично викону- 
ється на підвищеній розрядності, у відпо- 
відності зі значенням оцінки числа обумов- 
леності матриці, а також обчислюються 
оцінки достовірності отриманого розв'язку 
по). 

Якщо матриця не кваліфікується 
за(6) як  машинно  вироджена, але 
єдхсопдд 21, то СЛАР з наближеними 


вихідними даними, яку введено в комп'ю- 
тер, вважається некоректно поставленою 
при даній точності елементів матриці 
(інакше кажучи - виродженою в межах 
точності задання елементів матриці). До- 
стовірність обчисленого розв'язку не мож- 
на гарантувати. 

У цьому випадку користувач отри- 
мує комп'ютерний результат розв'язуван- 
ня задачі з оцінками його достовірності на 
основній розрядності та відповідне засте- 
реження. Крім того, розв'язування СЛАР 
продовжується на підвищеній розряднос- 
ті, у відповідності зі значенням числа обу- 
мовленості матриці. Обчислюються оцін- 
ки достовірності комп'ютерного розв'яз- 
ку. Збільшуючи довжину машинного сло- 
ва, у відповідності до оцінки числа обу- 
мовленості матриці, завжди можна отри- 


мати комп'ютерний розв'язок, який буде 
досить близьким до математичного роз- 
в'язку задачі. 

Для реалізації обчислень з довіль- 
ною розрядністю використано бібліотеку 
програм СМР (ОСХО Миїйріе-Ргесізіоп 
ІЛЬгагу) (21, 22|. Нижче наведено приклад 
дослідження СЛАР в Іпраг50їЇуег. 

Задача. Дослідити та розв'язати за 
допомогою Праг50іумег на СКІТ-4 систему 
лінійних алгебраїчних рівнянь Ах - б, де 


Аз ад)» ії -Ізп, поеЗучзі, 
у «1,2,... 
аз с по1-тахі(, )). 


ді --Нй-і 5; 
Отже, матриця системи має вигляд 
п-Іі п-1 п-2 .. 2 1 


п-Іі п-2 п-2 .. 2 1 
п-2 п-2 п-3 .. 2 1 

жа 
2 2 2 зо 
1 1 1 .. 1 0 


Елементи правої частини системи обчис- 
люються за формулами: 


раду, Р еп-і, якщо іх22; 
Брепоі-і, о якщої» 2. 


Порядок системи п « 1000. 
Точний розв'язок системи: 


х-(0 1 0 .. 0) 


Незважаючи на те, що елементи матриці 
та правих частин задано точно, при вве- 
денні їх в комп'ютер здійснюється заок- 
руглення до машинної точності тасПпер5 
ПО). Це враховується в процесі до- 
слідження таких задач в Іпраг5оїуег. 

Фрагмент протоколу розв'язування 
задачі: 


РЕОВІЕМ: зоіуїпд ої Рпе 1іпеаг аїдергаїіс 
зувкет мікр а зуттебіс 
робібіуе дебіпед паїргіх 


Ргосе55 ої іпуе5рідагіпд апа з5зоіуїпд 
ройріе ргесізвіоп 

МЕТНОГ: Сроїеї5ку десотровініоп 
ВЕ50О1Т58: 


ТНЕ МАТБІХ І5 МОТ РОЗІТІУЕ РЕКІМЕР !!! 


Митрег ої ргосез55ез: 1 
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МЕТНОГ: бацєз8 десотровіїіоп 


БЕЗОГТ58: 


ТНЕ МАТБІХ І5 МАСНІМЕ-5ІМСОГАК 1!!! 


РЕБЕСІЗІОМ: 128 


5ОБОТІОМ 
Гіге5і 4 соптропепі5 ої зоіціїіоп аге: 
0 1 0 0 


Сотригкасіопа) еггог іп Епе зоіціїіоп: 
0.00000еч00 


Тпе уесгог ої з5о01цііоп аге зиссез5іїціїу 
зсогед іп СБе Біїе гезців.опі 
Митрег ої ргосез85ез: 2 

З протоколу видно, що, оскільки 
матриця системи є симетричною, то проб- 
ним алгоритмом для дослідження задачі в 
Іпраг5оїмег | вибрано | алгоритм | І1/. 
розвинення як найбільш економічний за 
використанням обчислювальних ресурсів 
та часу виконання для таких матриць. Од- 
нак, в процесі дослідження цим алгорит- 
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мом на архітектурі 1 СРО та 1 СРО матри- 
ця виявилася недодатньо означеною й по- 
дальше дослідження задачі в Шраг50Їуег 
автоматично продовжено алгоритмом 1(- 
розвинення. 

У процесі дослідження СЛАР алго- 
ритмом 1У-розвинення, матриця виявила- 
ся машинновиродженою. Подальше роз- 
в'язування цим гібридним алгоритмом ав- 
томатично було продовжено на розряд- 
ності 128. Отримано очікуваний розв'язок 
задачі, видано 4 контрольних його зна- 
чень. Причому, для ефективного  до- 
слідження СЛАР цим алгоритмом автома- 
тично було вибрано топологію гібридного 
комп'ютера 2 СРО - 2 СРО. Весь вектор 
розв'язку збережено у файлі гезції. оп. 

Структура Шшраг5оЇуег 

На рис. | представлено структуру 
Іраг5оїЇуег: 


Користувач 


Інтелектуальний 
інтерфейс 


База даних 
та блок 
пояснень 


Алгоритмічно програмні База 
засоби дослідження та знань 
розв'язування задач 


Рис. 1. Структура ШрагЗоЇуег 


Складові частини Шраг5оЇуег: 
е інтелектуальний інтерфейс для спілку- 
вання з користувачами на мові пред- 
метної області; 


е- база знань з предметних областей, що 
розглядаються, на основі яких побудо- 
вано їх формальні моделі; 
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е база даних, до якої входять введені дані 
та результати досліджень задач в ком- 
п'ютері, що використовуються для 
адаптивного створення гібридного ал- 
горитму та синтезу відповідної програ- 
ми розв'язування; 

е блок пояснень, в якому накопичується 
інформація про хід обчислювального 
процесу задачі та результати аналізу 
достовірності отриманих розв'язків; 

е алгоритмічно-програмні засоби, за до- 
помогою яких реалізується автоматиза- 
ція всього процесу дослідження та роз- 
в'язування задачі на комп'ютері при 
використанні підвищеної комп'ютерної 
розрядності для забезпечення достовір- 
ності розв'язку задачі, а також засобів 
штучного інтелекту (нейронні мережі) 
для автоматичного розпізнання типу 
розріджених матриць. 

З метою побудови бази знань для 
кожної предметної області у вигляді від- 
повідних формальних моделей було про- 
ведено дослідження, що включає: 

е визначення типів задач з урахуванням 
різних обчислювальних схем, різних 
структур даних та їх властивостей; 

е. встановлення параметрів для кожної за- 
дачі, достатніх для реалізації автома- 
тичного дослідження її математичних 
властивостей та розв'язування; 

е аналіз складу необхідних алгоритмів та 
програм для розв'язування визначених за- 
дач з урахуванням архітектурних та тех- 
нологічних особливостей комп'ютера; 

е визначення логічних зв'язків та перехо- 
дів з метою реалізації автоматичного 
дослідження та розв'язування задач з 
наближеними даними; 

е аналіз системних компонентів (копілято- 
рів, процесорів СРО та СРО, оператив- 
ної пам'яті тощо), необхідних для орга- 
нізації обчислювального процесу на 
гібридному комп'ютері. 

Розроблено інтелектуальний інтер- 
фейс для спілкування користувачів з 
ІСКМ, який задовольняє таким вимогам: 


е інтерактивна підтримка на всіх етапах 
розв'язування задачі без додаткової 
документації; 

е врахування різних рівнів підготовки ко- 
ристувачів та їх завдань - передбачити 
можливість використання в автоматич- 
ному режимі дослідження та розв'язу- 
вання задачі (без втручання користува- 
ча) та в інтерактивному режимі, нада- 
ючи можливість користувачам прийма- 
ти участь у процесах дослідження зада- 
чі, вибору кількості процесорів тощо; 

е постійний контроль та перевірка несу- 
перечливості прийнятих рішень 
користувачем; 

е забезпечення розв'язання задачі від по- 
чатку до кінця або пояснення причин не- 
можливості обчислювального процесу. 

Діалог з користувачем забезпечує не 
тільки зручні способи задання даних задач, 
їх розв'язування та візуалізацію результа- 
тів задачі в Шраг5оіуег, але також сприяє 
поповненню знань з обчислювальної мате- 
матики, паралельних обчислень на гібрид- 
них комп'ютерах. Із цією метою надається, 
за бажанням користувача, детальна інфор- 
мація про технологію дослідження задач з 
наближеними даними та аналізу достовір- 
ності результатів, глосарій термінів з об- 
числювальної математики, низку прикладів 
задання вхідних даних. 

На програмному рівні Іпраг50імег - 
це сукупність програмних засобів взаємо- 
дії користувача з комп'ютером, що вклю- 
чає діалогові засоби, блок планування та 
управління обчислювальним процесом, 
програмні компоненти для вводу даних за- 
дач та виводу результатів тощо. 

Використання Праг5оЇуег 

Розроблена інтелектуальна система 
Іпраг5оімег функціонує на гібридному 
комп'ютері СКІТ-4 та може використову- 
ватися для математичного моделювання 
фізико-технічних процесів локально, в Ін- 
тернет або в Українській академічній грід- 
системі. 
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Так, система Іпраг50іуег та її окремі 
програми, що реалізують конкретні гібрид- 
ні алгоритми, були використані при мате- 
матичному моделюванні процесів в'язкого 
руйнування товстостінних елементів тру- 
бопроводів з дефектами зношення (спільно 
З Інститутом електрозварювання 
їм. Є. О. Патона НАН України) |23|. Резуль- 
тати математичного моделювання надають 
можливість визначати залишковий ресурс 
відповідальних конструкцій та приймати 
обгрунтовані рішення щодо подовження 
нормативних термінів їх безпечної експлу- 
атації. Зазначений методологічний підхід до 
чисельного моделювання не має світових 
аналогів і використовується в Інституті 
електрозварювання їм. Є.О. Патона НАН 
України для потреб атомної енергетики, а 
також трубопровідного транспорту 
України. 

Також алгоритмічно-програмні засоби 
Шраг5оіуег було використано в математич- 
ному моделюванні задачі стійкості нових 
композитних матеріалів із застосуванням 
тривимірної моделі «волокон кінцевих роз- 
мірів» для Інституту механіки 
імені С. П. Тимошенка НАН України |24- 
26). Задача зводиться до розв'язування част- 
кової узагальненої алгебраїчної проблеми 
власних значень для розріджених матриць, 
що була розв'язана новими гібридними 
алгоритмами методу ітерацій на підпросторі 
124,25). Обчислені мінімальні власні зна- 
чення дають можливість визначити вели- 
чини критичних параметрів стійкості 
шарувато-композитного матеріалу | при 
стисненні поверхневим навантаженням. 

У результаті розв'язування вказаних 
задач за допомогою створених високопро- 
дуктивних алгоритмічно-програмних засо- 
бів, отримано прискорення до 60-30 раз у 
порівнянні з результатами розв'язання 
традиційними програмними засобами на 
однопроцесорних комп'ютерах. 

Висновки 

Запропоновано інтелектуальну сис- 
тему комп'ютерної математики для до- 
слідження та розв'язування задач основ- 
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них класів обчислювальної математики з 

такими функціональними можливостями: 

е постановка задач обчислювальної мате- 
матики на мові предметної області; 

е звільнення користувачів від роботи по 
дослідженню математичних властивос- 
тей задач з наближеними даними, ство- 
ренню відповідних гібридних алгорит- 
мів та програм, аналізу достовірності 
результатів, що суттєво підвищує ефек- 
тивність математичного моделювання 
фізико-технічних процесів; 

е інноваційні способи ідентифікації порт- 
рету розріджених матриць на основі 
нейромереж, регуляризації різних струк- 
тур розріджених даних до стандартних 
видів для ефективного розв'язування 
задач; 

е. підвищення точності комп'ютерних роз- 
в'язків задач з наближеними даними, 
використовуючи багаторозрядну ариф- 
метику; 

е автоматизація всіх процесів розпарале- 
лення на гібридному комп'ютері при 
ефективному використанні комп'ютер- 
них ресурсів; 

е візуалізація та пояснення результатів 
розв'язування. 

Застосування інтелектуальної систе- 
ми Іпраг5оіуег та її складових у матема- 
тичному моделюванні прикладних додат- 
ків, що зводяться до розв'язування задач 
обчислювальної математики на гібридних 
комп'ютерах, дає можливість істотно по- 
кращити результати моделювання, скоро- 
тити час на його реалізацію 1, відповідно, 
заощадити кошти на проведення високо- 
вартісних натурних експериментів. 

Забезпечується розв'язування прик- 
ладних задач надвеликої розмірності (до 
десятків млн. ступенів свободи) і отриман- 
ня чисельних результатів з більшою до- 
стовірністю, завдяки застосуванню підви- 
щеної розрядності для погано обумовле- 
них та вироджених задач, в межах ком- 
п'ютерної розрядності. 
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